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ABSTRACT 

Navier-Stokes computational fluid dynamics calculations are presented for isolated, half-span, and full-span V-22 
tiltrotor hover configurations. These computational results extend the validity of CFD hover methodology beyond 
conventional rotorcraft applications to tiltrotor configurations. Computed steady-state, isolated rotor performance agrees 
well with experimental measurements, showing little sensitivity to grid resolution. However, blade-vortex interaction 
flowfield details are sensitive to numerical dissipation and are more difficult to model accurately. Time-dependent, 
dynamic, half- and full-span installed configurations show sensitivities in performance to the tiltrotor fountain flow. As 
such, the full-span configuration exhibits higher rotor performance and lower airframe download than the half-span 
configuration, Half-span rotor installation trends match available half-span data, and airframe downloads are reasonably 
well predicted. Overall, the CFD solutions provide a wealth of flowfield details that can be used to analyze and improve 
tiltrotor aerodynamic performance. 


NOTATION 

a speed of sound 

A rotor disk area, 7rR J 

c local chord length 

Cq rotor torque coefficient, Q//>(QR'f RA 

Ct rotor thrust coefficient, T/p(QR) 2 A 

DL/T airframe download divided by total thrust 
FM rotor figure of merit, C T -JC T /2/C q 
M Mach number, v/a 

M‘c„ blade section normal force coefficient times Mach 
number squared, N / -k pa 2 c 
N blade section normal force 

Q rotor torque 

r radial coordinate 

R blade radius 

Re Reynolds number at the rotor tip, p(QR)c tip /p 
T rotor thrust 

v local velocity 

x streamwise coordinate (+aft) 

y spanwise coordinate (+right) 

z normal coordinate (+up) 

F circulation 

0 blade collective angle at r/R = 0.75, degrees 

p air density 

G rotor solidity, NbC rcf /7tR 
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'F azimuthal angle, degrees 

CO vorticity, 1/sec 

Q rotor rotational speed, rad/sec 

INTRODUCTION 

Tiltrotor aircraft are recognized for their ability to 
significantly change both the military and civilian aviation 
transportation landscapes. The range and speed of a 
turboprop airplane is augmented by the ability to operate in 
and out of confined areas like a helicopter. For civilian 
operations this means reduced impact on an already 
overloaded airspace system and reduced infrastructure costs. 
For military operations, increased payload and range with 
reduced aerial refueling operations are possible when 
compared with helicopters currently performing the same 
remote area missions. The V-22 Osprey is the first 
production military tiltrotor aircraft. 

As with most new aircraft configurations, analysis tools 
need to be validated for regions beyond their conventional 
operation. In hover, tiltrotors differ significantly from 
helicopter rotors, which operate solely in edgewise flight, 
due to their highly twisted, low aspect ratio blades and 
high disc loading. These design aspects anse because the 
blades must operate in both propeller and helicopter rotor 
mode. Thick inboard airfoil sections display stall delay due 
to three-dimensional and centrifugal effects. In addition, 
airframe download prediction is more critical for tiltrotors 
than for conventional rotorcraft. The ability to accurately 
analyze and understand these features is a critical 
requirement for optimum tiltrotor design. For example, a 
0.01 change in figure of merit (FM), a measure of hover 
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efficiency, is equivalent to approximately 380 lbs. in useful 
load at constant power for a twin engine V-22 (C T = 0.014). 
In addition, about 1 0% of the rotor thrust is lost to airframe 
download for tiltrotor configurations. 

This study aims to investigate high-fidelity 
computational fluid dynamics (CFD) methods for predicting 
the performance of tiltrotors in hover, both for isolated 
rotors and full configuration aircraft. CFD methodologies 
have been investigated and applied extensively to 
conventional helicopter rotors. In this paper, they will be 
applied to two model rotor experimental test cases and to a 
full aircraft V-22 configuration. The aim is to extend the 
validity of the analysis methods, enabling investigation cf 
current full-scale tiltrotor aircraft, exploration of new design 
concepts, and development and better understanding of 
wind tunnel tests. 

Previous computational work on helicopter hover is 
extensive. Recent works in the field include Wake and 
Baeder [1], Beaumier, Pahlke, and Celli [2], and Pomin 
and Wagner [3]. Using the Navier-Stokes overset grid 
methodologies employed here, Strawn [4] showed good 
performance prediction for a 4-bladed UH-60 rotor. 
However, even with 64 million grid points for a quarter 
domain, the detailed wake geometry was not captured 
correctly, and grid convergence was only shown for 
integrated quantities. 

Previous work on tiltrotor configurations in hover is 
more limited, in part, due to a smaller number of tiltrotor 
configurations and experimental databases. An early Navier- 
Stokes simulation using a single block grid (240,000 grid 
points), momentum-source actuator disk, and a simplified 
half-span wing geometry [5] produced reasonable rotor 
performance, but download was significantly overpredicted. 
Numerous tiltrotor hover computations have been performed 
by Rajagopalan [6,7]. In these works, non-body-fitted 
Cartesian grids and time-averaged, momentum-source 
modeling of the rotor disk are used in efficient Navier- 
Stokes analyses of isolated (400,000 grid points) and 
installed (700,000 grid points) tiltrotor configurations. 
Isolated rotor performance is accurately predicted except at 
high thrust conditions, as are the incremental installation 
effects of half-span models. Meakin [8] computed the flow 
for a V-22 rotor in the presence of a half-span wing using 
the Navier-Stokes moving body, overset grid (2.5 million 
grid points) methodology employed here. Download was 
accurately predicted, but figure of merit was underestimated, 
although the comparison was made for slightly different 
rotors. 

Results in this paper will show validation of CFD for 
isolated tiltrotor hover performance by comparison with 
experimental test data. Effects of flow solver operation and 
grids will be investigated in detail to determine appropriate 
operational guidelines. Details of the three-dimensional 
isolated hover flowfield will also be studied - data not 


available from lower fidelity comprehensive codes or 
momentum-source models. The CFD methodology will 
then be applied to investigate the hover aerodynamics of 
half- and full-span tiltrotor aircraft configurations, including 
installed rotor performance, airframe download, and detailed 
flow physics. 

METHODOLOGY 

CFD calculations use the Reynolds-averaged Navier- 
Stokes computational fluid dynamics code OVERFLOW-D 
[9]. It is based on the OVERFLOW code, which was 
developed at NASA and has been applied to a wide range of 
fluid dynamics problems. OVERPLOW-D includes major 
modifications for rigid body motion of components. 
Solutions are computed on structured, overset grids using 
body-conforming “near-body” grids and automatically 
generated Cartesian “off-body” grids [10] in the wake and 
farfield. Several modifications to the code for calculations cf 
hovering rotors have been made [11], In particular, isolated 
rotor flows can be computed as steady-state problems in a 
blade-fixed reference Same with the addition of a rotational 
source term to the finite difference equations. 

More complex hover configurations require the time- 
accurate modeling of the moving rotor blades. User-defined 
subroutines prescribe the arbitrary six degree-of-freedom 
blade motion. Grid motion requires recalculation of the 
domain connectivity, including hole cuts and interpolation 
coefficients, at each time step, as the near-body grids move 
through the stationary off-body grids. Reuse of information 
from the previous time step enables order of magnitude 
speed-ups compared to domain connectivity solutions from 
scratch. Using this technique, the domain connectivity 
work can be efficiently performed in less than 20% of the 
time required for the flow solver. 

Solutions are computed on large parallel computers or 
a network of PCs/workstations communicating with the 
Message Passing Interface (MPI) protocol. Both the domain 
connectivity and flow solver modules have been parallelized 
for efficient, scalable computations using MPI [12]. Coarse 
grain parallelization on large numbers of processors is 
achieved by distributing grids among the processors, and, if 
necessary, splitting them as appropriate into smaller blocks 
to prevent bottlenecks. Boundaries that are created in the 
splitting process have explicit boundary conditions, similar 
to intergrid boundaries of the original grid system. 

ISOLATED ROTOR MODELING 
Isolated Rotor Test Cases 

A first step in the application of CFD for tiltrotor 
configurations in hover is the investigation of isolated 
tiltrotor Navier-Stokes CFD analysis capabilities. 
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Computational simulations were run for two experimentally 
tested isolated rotors, TRAM and JVX. 

The Tilt Rotor Aeroacoustics Model (TRAM) is an 
extensive wind tunnel model constructed to facilitate future 
tiltrotor aeromechanics research. The geometry is a 
0.25-scale V-22 right-hand nacelle and 3-bladed rotor with 
geometric and dynamic scaling. The isolated TRAM rotor 
was tested in the Duits-Nederlandse Windtunnel Large 
Low-speed Facility (DNW-LLF) in the spring of 1998. It 
provided a significant new source of aeroacoustics, 
performance, and structural loads data for validation of 
tiltrotor analyses. Among the aerodynamics data acquired 
were rotor performance and blade pressures. Hover runs were 
performed at nominal tip Mach numbers of 0.58 and 0.62. 
Rotor operational limitations resulted in the 0.62 hover 
data rather than the V-22 hover tip Mach number of 0.72. 
Hover testing was performed with the rotor shaft axis at 0 
degrees (helicopter mode) and -76 degrees (airplane mode). 
Support blockage and interference are reduced in airplane 
mode. Performance data in this mode are considered to be 
more accurate based on discussions by Johnson [13] and are 
used for comparisons with the computational results. 
Details of the test and data reduction are described in 
References 14 and 15. 

The JVX is a 0.658-scale model of an early V-22 
configuration with higher solidity obtained by increasing 
the chord 8.4% along the span. The isolated rotor was 
tested in the NASA Ames Research Center Outdoor 
Aerodynamic Research Facility (OARF) in 1984 [16]. 
Good data repeatability and reduced interference effects of the 
outdoor testing apparatus make this a high-quality dataset 
for computational validation. Data were obtained closer to 
the actual V-22 hover tip Mach number at 0.68 and larger 
tip Reynolds number of 5.9 million. 

The main physical characteristics of the two rotors are 
presented in Table 1, and nominal test parameters in hover 
are summarized in Table 2. Both rotors use V-22 airfoil 
sections and twist and thickness distributions along most cf 
the span (r/R > 0.25). Blade root fairings were slightly 
modified for each rotor. The large amount of nonlinear twist 
is typical of tiltrotors but significantly different from 
helicopter configurations. Details of the models and 
geometries have been described in Reference [13], 

Isolated Hover Analysis 

Rotor performance comparisons were made for both test 
cases. In addition, the TRAM rotor was used to 
investigate, in detail, the effects of grid topology, grid 
resolution, and flow solver algorithms as well as flowfield 
details. 

The 3-bladed TRAM rotor is shown in Figure 1. The 
inboard upper surface cuff geometry (r/R < 0.27) was 
obtained from a digitized scan of one of the model blades. 


Table 1. TRAM and JVX rotor physical characteristics 



TRAM 

JVX 

description 

0.25-scale V-22 
nacelle and rotor 

0.658-scale V-22 
development rotor 

rotor radius 

57 in. 

150.0 in. 

solidity, a 

0.105 

0.1138 

tip chord, c„ D 

5.5 in. 

15.79 in. 

airfoil 

V-22 XN series 

sections 

modified root 
fairing 

8.4% chord 
increase 

twist 

32 to -6 degrees, nonlinear 


Table 2. TRAM DNW and JVX OARF nominal test 


parameters in hover. 


TRAM 

JVX 

tip Mach number 

0.58, 0.62 

0.68 

tip Re (million) 

2.1 

5.9 

Ct/ct range 

0.05 - 0.17 

0.0-0.16 

collective range 

3-17 deg 

-4-17 deg 



Figure 1. Isolated TRAM geometry and baseline surface 
grids (every other point). 


while the remainder of the model used the CAD definition. 
The JVX is defined by CAD data. Elastic blade effects are 
minor in hover [13] and are not modeled. The centerbody is 
a V-22 spinner with a faired boattail. The hub is modeled 
without connecting flexbeam. The effect of this geometry 
simplification on the overall solution should be minimal 
due to low velocities in this region. The blade extends from 
0.I0R to 1.0R and is physically separated from the spinner 
by a distance of 0.026R. 
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Table 3. Grid Parameters. 



Baseline 

Alternative 

blade topology 

C-mesh: 

O-mesh: 


181x94x65“ blade 

121x94x65 blade 


1 1 1x67x65 hub 

81x67x65 hub 


109x47x65 tip 

88x49x65 tip 

near-body 

chordwise: 

2/3 of baseline: 

spacing 

.004-. 00 12c - LE 

271x139x97 blade 


0.0008c - TE 

166x101x97 hub 


normal: 0.0002 
spanwise: 
0.004R - root 
0.0024R - tip 
0.019R - max 

163x71x97 tip 

near-body outer 
boundary 


Ctip 

off-body 

spacing 

0. 10c tip 

0.05, 0.14c tlp 

outer boundary 

5R 


a chordwise x spanwise x normal 


Baseline Grid Generation 

The baseline grid systems for the TRAM and JVX 
rotors are similar, so only details of the TRAM 
configuration will be provided. In the Chimera 
methodology, overset, structured near-body grids are 
generated about the geometry. They extend approximately 
one tip chord (c,j p ) away from the body and include 
sufficient resolution to capture boundary layer viscous 
effects. The baseline surface grids are also shown in Figure 
1. They use C-mesh topology blades and tip caps. Grid 
spacing parameters in all directions are detailed in Table 3. 
The first four points away from the blade surfaces have a 
constant spacing, verified to produce a y+ < 1. These 
spacing parameters closely match those determined for 
accurate drag prediction on transonic fixed wings using 
OVERFLOW [17], 

Off-body Cartesian grid generation is automatically 
performed by OVERFLOW-D. The finest off-body spacing 
for the baseline grid is 0.10c tip . This level-1 grid surrounds 
the blades and extends ± 1.23R in x and y and -0.58, +0.23 
in z. It is manually specified in order to contain the entire 
wake. A total of five progressively coarser levels are 
generated out to the farfield boundary, which is placed at 5R 
in all directions from the center of the domain. The grid 
spacings differ by a factor of two between each mesh level. 
The baseline grid contains 15.9 million (M) points: 6.2M 
near-body and 9.7M off-body. Where grid points fall inside 
the geometry, hole cutting is employed to blank out these 
points. A cut through the grid system in Figure 2 shows 
the near- and off-body grids, hole cuts, and overlap. 



Figure 2. Slice through isolated TRAM volume grids 
(every other point), black - near-body, red - finest level 
off-body, blue - coarser level off-body 


Parametric Grid Variation 

In order to determine the sensitivity of calculated 
tiltrotor performance to grid effects, several alternative grids 
were generated for the TRAM configuration. Total numbers 
of points range from 11.1 to 37.4 million. For the blades, 
an O-mesh topology was generated with spacings similar to 
the C-mesh. If not for the O-mesh topology requirement to 
place five points on the blunt trailing edge, the C-mesh 
chordwise trailing edge spacings could have been relaxed 
slightly. Because of the reduced wake resolution, the 
O-mesh topology blades and caps have 26% fewer points. 
Additionally, O-meshes are significantly easier to generate 
so that their accuracy is worth investigating. 

Grid density was varied in both the near-body (NB) 
and off-body (OB) meshes. Coarser off-body grids with 
0.l4c t i P spacing in the finest level and correspondingly 
coarser spacings in the farfield were generated. A finer off 
body grid with 0.05c,j P spacing was manually generated 
inside the finest off-body grid level of the baseline grid. It 
was generated in cylindrical rather than Cartesian 
coordinates with extent 1.1 2R in r and -0.1 9R, +0.1 5R in 
z. Cylindrical coordinates were used so that the fine spacing 
would always be normal to the rotational flow. Azimuthal 
spacing is 0.37 degrees (0.067c tip at the blade tip). Finally, 
near-body blades and caps were generated with 1.5 times 
the number of points in all directions (2/3 of the baseline 
spacings) and combined with the fine cylindrical off-body 
grid. The different grid systems investigated are detailed in 
Table 4. 




Table 4. Grid systems. 


ISOLATED HOVER RESULTS AND DISCUSSION 


grid 

near- 

body 

spacing 

finest 

off-body 

spacing 

off-body 

points 

% 

# points 
(million) 

C-mesh 

baseline 

0.10c,,,, 

61 

15.9 

coarse off-body 

baseline 

0.14cti„ 

44 

11.1 

fine off-body 

baseline 

0.05c tip 

cylindrical 

74 

24.3 

fine near-/off- 
body 

2/3 

baseline 

0.05c, i P 
cylindrical 

47 

37.4 

O-mesh 

baseline 

0.10c, iD 

68 

14.2 


JVX Rotor Performance 

CFD solutions for the JVX rotor were run for a number 
of collective angles and a tip Mach number of 0.68 using 
the baseline grid system. Figures 3 and 4 compare 
calculated and experimental thrust, torque, and FM values. 
Excellent agreement is seen across the thrust range. The 
computed figures of merit generally underpredict the 
experimental values by less than 0.02 at all thrust levels. 


Flow Solver Details 


For spatial discretization, the baseline operation of 
OVERFLOW-D used 4 ,h -order central difference for the 
inviscid terms, 2 nt, -order central difference for the viscous 
terms, 4 lh -order artificial dissipation, no 2 nd -order artificial 
dissipation, and 2 nd -order grid metrics. All grids overlap 
sufficiently to allow double fringing, where two boundary 
points are interpolated from adjacent overlapping grids. 
This allows gradient as well as solution information to be 
transferred more smoothly and accurately between grids. 
The isolated hover solutions were run using the steady-state 
formulation. A characteristic boundary condition was used 
at the outer boundary where freestream values are 
approximated using a source-sink distribution [11], An 
LU-SGS time advancement algorithm is combined with 
local time stepping to speed convergence to a steady state. 

The Baldwin-Barth one-equation turbulence model is 
used in the near-body grids with thin-layer viscous terms 
added only in the computational direction normal to the 
surface. Turning on the thin-layer viscous terms in all three 
computational directions in the near-body grids produced 
no discemable change in the solution. Except as noted, the 
near-body domain is assumed fully turbulent with no 
transition region specified. Off-body grids as well as the 
centerbody were modeled using the inviscid. Euler equations 
in order to reduce the non-physical numerical dissipation of 
the wake. 

In addition to grid variation, flow solver parameter 
variations were investigated. The first involved changes in 
the 4 th -order dissipation coefficient. Second, for comparison 
with more complex applications, a moving grid, unsteady 
hover simulation was performed. This dynamic simulation 
was compared with results using the steady-state 
formulation. The time-accurate analysis used an implicit 
1 “-order Pulliam-Chaussee pentadiagonal algorithm in the 
near-body grids and an explicit 3 rd -order 3-stage Runge- 
Kutta scheme in the off-body grids. 



Figure 3. Isolated JVX rotor hover performance: 
torque vs. thrust. 



Figure 4. Isolated JVX rotor hover performance: 
figure of merit vs. thrust. 


685 



TRAM Rotor Performance 


CFD solutions for the TRAM rotor were run for a 
range of collective angles from 8 to 16 degrees with a tip 
Mach number of 0.625. Comparisons of computations with 
the TRAM DNW hover performance data are shown in 
Figures 5, 6, and 7, detailing torque, FM, and collective 
trends with thrust. The agreement for rotor performance is 
again quite good, however, the slope of the Ct-0 (Figure 7) 
curve is too large. There is a consistent overprediction of 
torque (~0. 000020) and underprediction of figure of merit 
(0.01-0.02) with the correct trends predicted across the 
collective range. 

The discrepancies are mostly due to the fully turbulent 
assumption. At these local chord Reynolds numbers, some 
laminar flow undoubtedly exists on the lower surface given 
the favorable pressures gradients there (see below for a 
discussion of blade sectional pressure distributions). Skin 
friction measurements on a full-scale XV- 15 rotor [18] 
indicate that there is a possibility for laminar flow on the 
upper surface as well at lower collective angles. The V-22 
and XV- 15 rotors have similar twist and thickness 
distributions, albeit different planforms (tapered vs. straight) 
and airfoil sections. The XV- 15 test tip Mach numbers of 
0.56 and 0.69 equate to higher tip Reynolds numbers of 
4. 0-5.0 million. At collectives angles between 7 and 14 
degrees, upper surface transition was seen to recede from 
30% to 3% chord on average along the span of the blade. 

Based on an average location of the start of adverse 
pressure gradients in the pressure distributions, a laminar 
region was specified for the 14-degree collective case. The 
region extends from the leading edge to x/c = 0.60 on the 
lower surface and to x/c = 0.05 on the upper surface. This 
transition case is plotted in Figures 5, 6, and 7 and detailed 
in Table 5. It produced a 0.014 increase in FM due in part 
to a 0.000016 decrease in the viscous component of the 
torque coefficient. In this case, the viscous torque 
contribution (.000074) is less than 5% of the total torque. 
As a result, accounting for boundary layer transition further 
improves the performance comparisons to within the range 
of experimental data uncertainty and should also improve 
the correlation at other collective angles. 

Blade Sectional Quantities 

Blade pressures from the TRAM test were obtained at 
seven radial stations during hover testing in helicopter 
mode and integrated to obtain section normal force 
distributions. Experimental and computed pressures at 
radial stations r/R = 0.82 and 0.90 are compared in Figure 
8. Agreement is reasonable, but at higher thrusts, upper 
surface pressures are consistently overpredicted. At this 
point, it should be noted that integration of the 
experimental pressures underestimates the balance thrust by 



Figure 5. Isolated TRAM rotor hover performance: 
torque vs. thrust. 



Figure 6. Isolated TRAM rotor hover performance: 
figure of merit vs. thrust. 



Figure 7. Isolated TRAM rotor hover performance: 
thrust vs. collective. 
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Table 5. TRAM isolated rotor performance comparison at 14 degrees collective. 


Performance Wake geometry (first BVI) 


grid 

I Ct 

Co 

FM 

r/R 

z/c 

radius/c 

1 03 |max 

T/QRc 

C-mesh (baseline) 

f .01495 

.001660 

.779 

.86 

-.50 

.85 

.017 

.49 

transition 

1 .01513 

.001660 

.793 r 






coarse off-body 

; .01483 

.001659 

.770 

.86 

-.55 

.94 

.011 

.48 

fine off-body 

.01495 

.001675 

.772 - 

.87 

-.53 

.55 

.037 


fine near-/off-body 

i .01498 

.001681 

.771 

.87 

-.53 

.52 

.040 

.44 

O-mesh 

1 .01408 

.001562 

.756 

.87 

-.53 

.82 

.017 

.45 

dynamic 

; .01491 

.001705 

.755 

.86 

-.59 

.83 

.016 

.46 


at least 15%. Errors in surface pressure data or fidelity of the 
surface geometry may be the cause. Experimental normal 
force distributions are compared with computation in Figure 
9, which serves to show that the trends with r/R and 
collective are representative. Unfortunately, it is not clear 
from the pressure and normal force coefficient comparisons 
whether the calculated spanload distribution outboard of 
r/R = 0.70 is correct. The loads on this part of the rotor are 
largely due to the first blade-vortex interaction (BVI). CFD 
calculations tend to overpredict the effects of the first BVI 
due to excessive dissipation of the preceding tip vortex. 
The influences of blade-vortex interaction and vortex 
dissipation on spanload are discussed below in relation to 
grid effects. 

Flowfield Details 

Detailed investigation of the TRAM CFD solutions 
indicates mild flow separation on the cuff and inboard part 
of the blade at all collective angles computed. The separated 
region, as determined by a contour of zero skin friction in 
the rotating reference frame, is shown in Figure 10 for 14 
degrees collective. Surface streamlines in this reference frame 
are also shown and indicate three-dimensional relief is 
provided by the centrifugal effects forcing the flow outboard. 
The normal extent of the reversed flow is quite small, less 
than 3 percent of the local boundary layer thickness. These 
observations qualitatively match skin friction measurements 
on a full-scale XV- 15 rotor at similar collective angles and 
Cj/'g [18]. Despite rotor geometric differences, separation 
locations at r/R = 0.17, 0.28, and 0.50 and the overall 
spanwise flow region are in generally good agreement. The 
cause of the separation is related to the trade-offs in tiltrotor 
design requiring high root twist when operating as a 
propeller compared with desired conventional helicopter 
root twist. Based on a spanwise torque distribution, there 
does not appear to be any performance penalty associated 
with the inboard separation. Further investigation of the 
inboard flowfield may be useful for correlating stall delay 
characteristics for comprehensive codes, which use modified 
airfoil tables, since this phenomenon is directly modeled 



x/c 



Figure 8. Isolated TRAM rotor chordwise pressure 
distributions, symbols - test, line - calculated 



r/R 

Figure 9. Isolated TRAM rotor normal force 
distributions, symbols - test, line - calculated 



Figure 10. Isolated TRAM rotor upper surface streamlines 
(red) and separated region (black), Ct = 0.015. 




coarse OB grid baseline grid fine NB/OB grid 

Figure 11. Isolated TRAM rotor wake vorticity at blade, blue - low, red - high 


by CFD. A small separation region is also seen outboard at 
the foot of the shock near the leading edge. 

Grid Effects 

Performance results on the alternative grids are also 
shown in Figures 5, 6, and 7 and tabulated in Table 5 for 
computations at 14 degrees collective. The O-mesh 
underpredicts C T and FM compared with all the other grids. 
Flow visualization of the O-mesh results (not shown) 
indicates an immediate degradation of the wake vortex sheet 
due to the poor resolution behind the trailing edge. 
Additionally, the trailing edge details are not captured as 
accurately by the O-mesh. Reduced pressure on the lower 
surface aft of x/c = 0.8 along the entire span explains the 
underprediction of Ct. The spanloads and wake 
visualizations for the 0- and C-meshes are the same, 
indicating that the dominant problem may be the poor 
trailing edge treatment rather than the loss of the wake sheet 
from the trailing edge. Despite their ease of generation, 
O-meshes are not recommended for accurate rotor blade 
calculations with OVERFLOW-D. 

For the C-mesh grid systems, the effects of grid density 
on overall performance are minimal. Between the baseline 
and fine near-body grids, the viscous forces are the same, 
indicating grid convergence in the normal distributions. 
Differences are more apparent in the spanloads and wake 


structure. Table 5 also shows details of the first BVI. This 
first vortex encounter has a major effect on the flowfield in 
the outer blade region. 

Figure 1 1 shows the wake vorticity contours for the 
three off-body grid densities. For the baseline grid, the first 
BVI vortex has a radius that is 63% larger, a peak vorticity 
that is 58% lower, and a miss distance that is 0.03c, ip 
higher compared with the fine grid solution. The total 
circulation is the same for the two grid densities, indicating 
vorticity conservation despite the numerical dissipation. 
The baseline grid maintains the tip vortices until a wake 
age of 360 degrees, at which point they move into the next 
level coarser off-body mesh. 

Figure 12 shows the effects of vortex' dissipation on the 
spanload. In the baseline solution, a more diffuse and closer 
vortex induces higher local flow angles outboard of the BVI 
location (r/R = 0.86) and lower local flow angles inboard. 
The result is a peakier spanwise loading .distribution 
compared to the fine grid. This is supported by 
examination of pressure distributions. The coarse grid result 
shows combined effects, since the vortex is even more 
diffuse, significantly weaker, but has a larger miss distance. 
These grid convergence results agree with previous findings 
for helicopter rotors in hover [4], Convergence of integrated 
quantities can be achieved with careful attention to grid 
details, however, convergence of flowfield details has not 
yet been realized. 
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Figure 12. Isolated TRAM rotor normal force 
distributions for different grid densities. 


Flow Solver Effects 

Variation of the 4 th -order dissipation coefficient showed 
that there is no change in the solution for values less than 
the baseline value used here of 0.02. However, larger values 
adversely affect the performance predictions. Note that a 
dissipation coefficient of 0.04 is the default for 
OVERFLOW. Smaller values tend to make the solution 
more unsteady, but not unstable. 

Steadv-State and Dynamic Formulations 

For comparison with complex, moving body problems 
for which there is no transformation to a steady-state 
reference frame, a dynamic, moving grid, hover simulation 
was run. Exactly the same grid system and hole cutting 
methodology are used. The major difference is that the near- 
body blade grids move dynamically and time accurately 
through the stationary off-body grids. The moving grid 
calculation uses 3200 iterations per revolution, 0.1125 
degrees per iteration. 

For complex configurations with multiple rotors, a 
source-sink characteristic boundary condition may not be a 
good model for the bottom outflow boundary. In this case, 
the outer boundary condition uses freestream (M„ = 0) 
characteristics. Solutions using source-sink and freestream 
characteristic boundary conditions showed negligible 
differences when the outer boundaries are at least 5R away. 

Performance comparisons (Table 5) show that the 
calculated thrust is the same between the steady and 
dynamic calculations, but the torque is larger for the 
dynamic formulation, resulting in a 0.024 reduction in FM. 
This is due to a difference in the miss distance of the first 
BVI and corresponding effects on the outboard blade 
loading. 


In theory, there should be no difference in the solutions 
based on the choice of reference frame for the problem. In 
practice, however, several subtle differences in the 
methodology can arise. First, the time advancement scheme 
adds additional dissipation and dispersion. A time step 
convergence study was performed, and AT 7 = 0.1125 deg 
was determined to be sufficient for temporal accuracy, in 
agreement with the findings of other researchers [3]. No 
change in performance was noted for A*F values as small as 
0.0281 deg. However, due to the interaction of numerical 
diffusion between the temporal and spatial schemes, it may 
be necessary to perform a convergence study simultaneously 
in time and space to state this conclusively. 

Changing the order of accuracy of the temporal scheme 
in the off-body grids from the baseline 3 rd -order explicit to 
l^-order implicit showed little difference in overall 
performance. Lagged boundary condition information also 
affects the temporal accuracy. Grid splitting for parallel 
processing introduces numerous explicit Chimera 

boundaries into both the near-body and off-body grids. 
While temporal accuracy may be a cause for the steady-state 
vs. time-accurate discrepancies, it was not possible to 
produce the required magnitude of changes with the options 
investigated. 

Second, the steady simulation uses an analytical 
representation to model the rotational motion. The grid 
speed calculation in the steady formulation is exact while 
the dynamic formulation uses backward differences to 
compute the grid velocities. The steady simulation also 
uses an analytical rotational source term. This centrifugal 
force term should balance the right hand side residual from 
the unsteady formulation. Due to numerical and 
discretization errors, the terms may be slightly out cf 
balance, perhaps manifesting itself in a lack of conservation 
of angular momentum. 

At this time, the discrepancies between steady-state and 
dynamic hover results cannot be explained. However, they 
have significant repercussions for the ability to analyze 
complex simulations and comparison back to isolated 
performance based on a steady formulation. Similar 
discrepancies were seen elsewhere [3] and remained 
unresolved. 

Computational Convergence 

Solutions were run using 48-64 processors of an IBM 
SP3. For steady-state simulations and depending on the 
configuration and collective, convergence of thrust and 
torque is reached after 15,000 to 30,000 iterations, requiring 
approximately 20 to 40 wallclock hours. The values are 
oscillatory within ±0.5%. The magnitude of the 
oscillations is reduced by turning off local time-stepping 
after initial convergence. The dynamic simulation includes 
a 16% penalty due to recalculation of the domain 
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connectivity every iteration. Approximately 10 rotor 
revolutions are required to obtain converged force/moment 
values, equating to 32,000 iterations. Sample convergence 
histories for isolated steady-state (restarted from a lower 
collective angle) and dynamic simulations are shown below 
in Figure 14. 

INSTALLED ROTOR APPLICATION 

Investigation of isolated, hovering tiltrotors has shown 
that accurate rotor performance predictions (Figures 4 and 6) 
can be obtained given the guidelines of grid resolution and 
flow solver operation. This suggests that the CFD methods 
can be successfully applied to simulate more complex 
tiltrotor configurations, and, in particular, to investigate 
rotor installation effects. 

Half-span JVX Test 

The JVX rotor was tested in the presence of a 
0.658-scale half-span V-22 wing, simulated fuselage fairing, 
and image plane. This occurred during the same test as the 
isolated configuration. Because there are no tunnel-to-tunnel 
installation effects, the configuration increments from this 
test should be of especially high quality. Wing pressures 
were recorded at six span stations and integrated to obtain 
download measurements. Comparisons with the TRAM 
rotor are made at consistent Cj/o. Details of the model, 
test, and results can be found in References 16 and 19. 
Nominal test parameters in hover are the same as isolated 
(Table 2). 

Installed Hover Analysis 

A CFD model has been constructed of an installed 
TRAM rotor on a V-22 airframe configuration in hover. 
The blade geometry is the same used for the isolated hover 
analysis. The aircraft fuselage, wing, nacelles, and tail are 
actual V-22 geometry. There are some differences between 
the JVX half-span model described above and the full-span 
V-22 as modeled here. The non-dimensional rotor 
separation distance (separation/2R) for the V-22 is 1.226- 
while for the JVX it is 1.210. The JVX fuselage is a cmde 
representation and nacelle/wing fairings are absent. 
Compared with the flight article, the CFD model has 
approximate nacelle/wing fairings, simplified flaperon cove 
seals, and no protuberances (VGs, strakes, probes). In the 
CFD model, the flap has been faired into the wing to 
completely seal the gap. The wing flaperon angle is set to 
67 degrees and compared with data at the same setting. For 
simplicity, the nacelle inlet and exhaust locations have been 
modeled as solid surfaces rather than flow-through. Both 
full-span and half-span with image plane time-dependent 
CFD simulations were ran. The intent is to investigate the 


differences between half- and full-span configurations, as 
might be tested in a wind tunnel. 

Grid Generation and Flow Solver 

The grids and flow solver operation are almost 
identical to that used for the dynamic isolated baseline 
analysis. Differences are noted here. The blade grid system 
used in the isolated hover analysis was installed directly on 
the V-22 with a nacelle angle of 90 degrees. Rotor rotation 
on each wing is such that the blades pass from wing leading 
to trailing edge. The distance between the spinner and blade 
root is smaller than for the isolated configuration. The 
nacelle spinner rotates with the blades, sliding over a 
portion of the nacelle spinner collar. 

For the installed configurations, all near-body grids 
including the nacelles and fuselage are modeled as viscous. 
Viscous grid generation for the fuselage components relaxed 
the normal grid spacing due to significantly lower flow 
velocities and reduced wall-bounded viscous effects. The 
outer boundary spacing of the fuselage near-body grids 
matches with the finest off-body grid spacing, which 
remained at O.lOcnp. The finest level off-body grids enclosed 
not only the rotor plane as in the isolated analysis but also 
the complete airframe. Figure 13 shows a streamwise cross- 
sectional cut through the half-span volume grid system. 
The full-span grid system is approximately symmetric. The 
total number of grid points in the full-span model is 47.6 
million with 63% in the off-body grids. The half-span 
model contains 23.6 million points. The 4 th -order spatial 
scheme is used throughout. Double fringing is used in the 
blade and off-body grids while the fuselage near-body grids 



Figure 13. Slice through half-span V-22 volume grids 
(every third point), black - near-body, red - finest level 
off-body, blue - coarser level off-body 
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Figure 14. Isolated and installed TRAM rotor performance 
time histories at 10 degrees collective. 


only have enough overlap for single fringing. Time 
advancement schemes and time step are unchanged. The 
installed simulations apply freestream characteristic 
conditions to the outer boundary. The half-span simulation 
has an inviscid plane of symmetry boundary condition at 
the centerline. 

Both simulations were run for at least 31 rotor 
revolutions on 124 processors of an SGI Origin 2000. Each 
blade revolution requires 14.7 (full-span) and 8.6 (half-span) 
wallclock hours for 3200 iterations per revolution. Domain 
connectivity accounts for 20% of the wallclock time. 

INSTALLED HOVER RESULTS AND DISCUSSION 

The collective angle for the computed installed cases is 
10 degrees with tip Mach number and Reynolds number 
unchanged at 0.625 and 2.1 million, respectively. Due to 
previously noted discrepancies in absolute performance 
between steady-state and dynamic isolated hover analyses, 
an isolated rotor dynamic simulation at 10 degrees 
collective was run for comparison with the installed 
configurations. Configuration increments from consistent, 
dynamic analyses should be unaffected. 

Rotor Performance 

Rotor performance for the isolated, half-span, and full- 
span configurations are shown in Table 6 along with 
experimental data from the isolated TRAM DNW and JVX 



Figure 15 Installed vs. isolated rotor figure of merit 
increments. 


Table 6. Performance comparisons of installed V-22 
configurations. Calculations at 10 degrees collective. 


configuration 

C T /a 

C 0 /a 

FM 

AFM a 

isolated - steady 

.1002 

.00967 

.751 


isolated - dynamic 

.1005 

.00999 

.730 


half-span 

.0978 

.01021 

.688 

-.038 

full-span 

.1000 

.01002 

.724 

-.006 

DNW TRAM 

.1005 

.009505 

.768 


JVX isolated 

.1005 

.009533 

.797 


JVX half-span 

.1005 

.009921 

.766 

-.031 


a installation AFM estimated at constant Ct 


OARF tests. All FM increments due to rotor installation 
effects are estimated at constant thrust. Computational 
values for the installed configurations are averaged starting 
from the 20 th rotor revolution. Time histories of thrust and 
torque are shown in Figure 14 for steady-state and dynamic 
isolated, half- span, and full-span configurations. Time- 
dependent solutions indicate a three per revolution 
oscillation with the largest oscillations for the half-span 
rotor. 

The half-span analysis shows a 0.00028 loss in Ct 
( 2.7%) at constant collective and a 0.038 reduction in FM 
compared with the isolated calculation. At the same Cr/a, 
the JVX data for a half-span model without fuselage 
compared to an isolated rotor indicates a 0.031 FM 
reduction. The full-span CFD analysis shows almost no 
loss in thrust (< 0.5%) at constant collective and a much 
smaller performance penalty of 0.006 in FM compared with 
the isolated calculation. Hover performance increments due 
to rotor installation are summarized in Figure 15. Overall, 
the CFD calculations seem capable of correlating half-span 
installation increments with the available data. Comparison 
with accurate full-span increments should prove interesting. 
Absolute performance values are consistently underpredicted 
due to unresolved issues in the time-accurate simulations, 
although correlation improves somewhat when boundary 
layer transition corrections are applied. 


691 


Hi 


IT 


\ \ yj 
VVvUUVt'i ! it* L 


„.,v^ 


Vi * VV i 1 i f 4 — l •:. . X\% k\ 1 L V* J-i it*,.? X * %\ l * f f f * S f J l '-i 

■ v^HUVvii:'. " --- ' ’ ■ ■ * , rn '•i" 1 "-^ — - • •^-^ 1 Iii{//P1VM ■;• 

/~ T "r — ^r-T— * 3-- -1* 





p: / - 


half-span 



Figure 16. Time-dependent particle traces through V-22 simulations, 
colored by particle release time, blue - earliest, red - latest. 


Flowfield Details 

It is well known that there is an installation penalty 
due to a fountain of recirculating flow formed by the 
presence of the wing and second rotor, represented as either 
an image plane or full-span configuration. Some of the rotor 
downwash over the wing is deflected spanwise and inboard 
with the majority spilling over the wing leading and 
trailing edges. When the spanwise flow hits the centerline, 
it is deflected upward due to the image plane or spanwise 
flow from the opposite wing. The flow is then reingested 
into the rotor, reducing its performance as it passes over the 
wing. 

Based on the time histories of the full-span analysis 
(Figure 14), the left and right rotor flowfields are initially 
symmetric and start becoming asymmetric after six rotor 


revolutions. Flow visualization indicates that the fountain 
is initially biased towards the left rotor but then switches to 
the right. The average calculated FM difference between the 
two rotors is only 0.002. The asymmetry has been noted 
both experimentally [20] and computationally [7] with a 
frequency 1-2 orders of magnitude smaller than the rotor 
rotation frequency, which is just beyond the length of the 
current simulation. 

The differences in half- and full-span fountain flows can 
be seen in the particle traces of Figure 16. Particles were 
released in a line connecting the rotor hubs and traced over 
three rotor revolutions. Due in part to the swirl component, 
the fountain is highly three-dimensional. The half-span 
calculation indicates a significantly larger fountain height 
and spanwise extent compared with the full-span 
calculation. The large spanwise extent is seen both above 
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the rotor plane and in the disturbed flow behind and below 
the wing. Away from the fuselage, the organized wake of an 
isolated rotor is seen. 

Representative instantaneous vorticity contours in a 
streamwise plane at the rotor centerline are shown in Figure 
17. In the fountain region the tip vortices are initially 
convected towards the wing and then deflected upwards in 
the fountain. Several tip vortices can be seen in the 
recirculating flow. In the full-span case, the right rotor 
ingests the tip vortices from both rotors at this time step. 
For the half-span case, the larger size of the fountain and 
more disorganized region around the wing are easily seen. 
These features equate to a larger performance penalty 
compared with the full-span simulation. The blade wake 
can also been seen convecting down onto the wing. Overall, 
the characteristics of the two fountain flows match the 
experimental flow visualizations by Polak [20]. 

The recirculating flow and wing presence is evidenced 
in the blade loads in Figure 18. Away from the fuselage at 
azimuthal angles ( l F) of 30 deg and 150 deg (0 deg along 
the streamwise axis and increasing in the direction of 
rotation), the loads are similar to the isolated analysis. 
There is no difference between the half- and full-span loads 
at these two locations. Over the fuselage ('P = 270 deg), 
large reductions in outboard loading, especially for the half- 
span case, indicate significant downwash from the 
recirculating region. This results in an overall thrust loss. 
Loadings inboard on the rotor blade while over the wing 
show a small increase due to the simulated sround effect. 


There is residual inboard loading remaining at 30 deg 
azimuth. In the full-span simulation, this may counteract 
the reduced tip loadings over the wing enough to explain 
the small installed average thrust loss. It does not 
compensate for the even larger outboard thrust loss over the 
wing in the half-span case. 

Rotor-on-Wing Effects 

Airframe downloads for the installed configurations are 
shown in Table 7 along with experimental data from the 
JVX OARF test ( Ct/g = 0.100). The calculated average 
download includes all surfaces except the nacelle inlet and 



Figure 18. Isolated and installed TRAM rotor normal 
force distributions. 



half-span full-span 

Figure 17. V-22 fountain vorticity contours, blue - low, red - high 
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exhaust planes, which while modeled as solid are more 
realistically set to freestream pressure and accounted for in 
jet thrust. Values are averaged starting from the 20 lb 
revolution. Time histories of airframe forces are shown in 
Figure 19. Calculations indicate a rearward force that is 
trimmed out in flight with forward stick and/or nacelle tilt. 
Average side force for the full-span configuration is 
essentially zero. 

The half-span CFD calculation shows increased 
download/thrust (0.130) compared with the full-span 
calculation (0.100). Download breakdown by component is 
also shown in Table 7, where the wing component extends 
only to the side of body (r/R < 1.01). The two CFD 
configurations carry the same download on the wings, so an 
additional load on the fuselage and sponsons accounts for 
the increased half-span total. The extra relief that is 
provided by the full span to reduce the fountain flow results 
in lower download, in addition to the previously discussed 
smaller rotor performance penalty. Taking into account 
configuration differences, the calculated wing component 
download (0.069) is in reasonable agreement with the JVX 
data which shows a wing download of 0.088 out of a total 
0.106 (±0.01). The experimental wing download is 
calculated from integrated wing pressures only up to the 
side of body. The total includes an extrapolation to the 
image plane, accounting for the fuselage. 

In Reference 21, the download breakdown between 
wing and fuselage on a semispan V-22 model with image 
plane was reported as 0.062 (62%) and 0.038 (38%), 
respectively, out of a total 0.100. This data indicates that 
the half-span fuselage contribution from CFD is 
overpredicted. It is theorized that a viscous image plane 
might reduce the high velocities in the fountain flow. This 
would have the effect of reducing the fuselage download as 
well as increasing rotor performance. Nevertheless, overall 
download trends are reasonably well predicted by the CFD 
calculations given the scatter in available experimental 
data. The higher download on the half-span model is 
consistent with the physics of the computed fountain flows. 

Figure 20 indicates the time-averaged pressure force in 
the download direction for the full-span V-22 calculation. A 
zero contour line separates regions of up- and download. 
The viscous contribution is negligible. The majority of the 
wing download occurs at the wing tips. A large upload 
region is seen at the flap hinge line and leading edge (not 
visible). Most of the fuselage download contribution is due 
to consistently negative pressures on the underside. Except 
in the fountain flow region, there is an upload on the 
topside of the fuselage. Flow spillage off the wing leading 
and trailing edges at the wing/fuselage fairing produces two 
download regions on the upper surface of the sponsons. 
Flow acceleration over the sponson sides creates an upload. 
The sponsons may be a worthwhile area of investigation for 
download reduction. 


Table 7. V-22 download comparisons and breakd 
by airframe component, Crfct = 0.100 (-. j s upl 0a j" n 


component 

I DL/T 

“ r 

! DL/T 

% : 

: DL/T V- 

wing 

! .088 

! .069 

53 


wing+ a 

: ,io6 i 

fuselage 


! .034 

26 

•° 15 ~TT 

sponsons i 


i .022 

17 


nacelles 1 


1 .004 

3 

•005 

tail i 


j .ooi 

1 

: -.002 “T- 

TOTAL i 

■i .130 i 

.100 

a . , - . , - — c 



Figure 19. V-22 airframe force time histories. 





Figure 20. V-22 time-averaged airframe download 
pressure component, red - download, blue - upload, 
black line - zero contour 
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CONCLUSIONS 


REFERENCES 


Navier-Stokes CFD has been applied to isolated, half- 
span, and full-span tiltrotor hover configurations, focusing 
on rotor performance, flowfieid details, and airframe 
download. Guidelines have been given for grid generation 
and flow solver operation. The following conclusions are 
made from the data presented: 

• Isolated tiltrotor performance across the thrust range is 
accurately predicted using a steady-state CFD formulation 
and taking into account boundary layer transition. 
Integrated forces are relatively insensitive to grid density, 
however, blade load distributions and wake geometry 
still show effects of numerical diffusion and cannot yet be 
validated. 

• Flowfieid and performance differences exist between 
steady-state and time-dependent, dynamic isolated hover 
solutions. The unresolved discrepancies have 
implications for the ability to calculate complex, installed 
tiltrotor configurations and comparison with isolated 
performance by steady-state calculations. Isolated rotor 
performance was consistently underpredicted by dynamic 
simulations, although incremental effects should not be 
affected. 

• Calculations compared rotor performance and airframe 

download between half- and full-span V-22 

configurations. Differences can be traced to variations in 
the fountain flows. Download calculations show 
reasonable agreement with experimental measurements. 
The half-span rotor performance installation increment 
matches available test data. 

• High-fidelity CFD methods provide a broad range of 
flowfieid details, such as rotor sectional data, separated 
regions, and download distributions, that can be 
investigated to obtain improved performance from 
tiltrotor designs. 

Future work is required and will focus on investigating 
the discrepancies discussed and expanding the CFD 
capabilities beyond performance prediction. 
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